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Abstract. We discuss the quantization of sound waves in a fluid with a linear dispersion 
relation and calculate the quantum density fluctuations of the fluid in several cases. These 
include a fluid in its ground state. In this case, we discuss the scattering cross section of light 
0^ , by the density fluctuations, and find that in many situations it is small compared to the thermal 

■ fluctuations, but not negligibly small and might be observable at room temperature. We also 
consider a fluid in a squeezed state of phonons and fluids containing boundaries. We suggest that 
the latter may be a useful analog model for better understanding boundary effects in quantum 
field theory. In all cases involving boundaries which we consider, the mean squared density 
fluctuations are reduced by the presence of the boundary. This implies a reduction in the light 
scattering cross section, which is potentially an observable effect. 

O 
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1. Introduction 

It is well known that quantized sound waves, whose excitations are phonons, share several 
properties with relativistic quantum fields, such as the electromagnetic field. Here we will be 
primarily concerned with fluctuations in the phonon ground state. These fluctuations can share 
certain features with electromagnetic vacuum fluctuations, including the creation of Casimir- 
type forces pQ. The forces due to quantum or stochastic sound fluctuations have been discussed 
recently by several authors [2j [31 [H El [TJ. In this paper, we will study quantum density 
fluctuations in a fluid. The fluid in question could be either a classical or a quantum fluid, 
but we will consider only the regime where the dispersion relation is approximately linear. Thus 
the analog models which we consider are rather different than those which model the Hawking 
process El IIP] , The topic of hydrodynamic fluctuations has been reviewed, for example, in 
Refs. [HJ [12]. In Sec. El we will discuss density fluctuations in a fluid in the phonon ground 
state without boundaries, and summarize work recently reported in Ref. [13] . It was argued 
there that light scattering by quantum density fluctuations might be observable, even in fluids 
at room temperature. The effects of a squeezed state of phonons will be briefly discussed in 
Sec. El In Sec. [U we turn to the issue of the effects of boundaries in the fluid, and argue that 
this can be a useful analog model for the effects of boundaries on relativistic quantum fields. 
We also suggest that the boundary effects, at least in principle, create observable effects in the 



fluid case. We report on the result of several specific calculations for different geometries. The 
results will be summarized and discussed in Sec. [5j 



2. Phonons and Density Fluctuations in a Fluid 

2.1. Quantization and the Density Correlation Function 

We consider the quantization of sound waves in a fluid with a linear dispersion relation, Vt q = cs q, 
where Q q is the phonon angular frequency, q is the magnitude of the wave vector, and cs is the 
speed of sound in the fluid. This should be a good approximation for wavelengths much longer 
than the interatomic separation. Let po be the mean mass density of the fluid. Then the 
variation in density around this mean value is represented by a quantum operator, p(x, t), which 
may be expanded in terms of phonon annihilation and creation operators as |14j 



P(x,t) = £(& q / q + 6t/*), (1) 
q 

where 

Here V is a quantization volume. The normalization factor in Eq. ([2]) can be fixed by requiring 
that the zero point energy of each mode be ^Ml q and using the expression for the energy density 
in a sound wave, 

U = ^p 2 . (3) 
Po 

In the limit in which V — > oo, we may write the density correlation function as 

(p(x, t) p(x', tf)> = / d\ Q g e^'Ax-O, At) f (4) 

where Ax = x — x' and At = t — t'. The integral may be evaluated to write the coordinate space 
correlation function as 

w n ~, , a> hpo Ax 2 + 3c%At 2 

WM ^, 0) ___^__^L_. (5) 

This is of the same form as the correlation function for the time derivative of a massless 
scalar field in relativistic quantum field theory, (</3(x, t) <£>(x', t')). (This analogy has been noted 
previously in the literature. See, for example, Ref. [33].) Apart from a factor of po, these two 
quantities may be obtained from one another by interchanging the speed of light c and the speed 
of sound cs- If c — > cs, then 

(V3(x, t) <p(x', t')) po (p(x, t) p(x', t')) . (6) 
In the limit of equal times, the density correlation function becomes 

Thus the density fluctuations increase as |Ax| decreases. Of course, the continuum description 
of the fluid and the linear dispersion relation both fail as |Ax| approaches the interatomic 
separation. Also note the minus sign in Eq. ([7]). This implies that density fluctuations at 
different locations at equal times are anticorrelated. By contrast, when cg|At| > |Ax|, then 



(p(x,t) p(x' ,t)) > and the fluctuations are positively correlated. This is complete analogy 
with the situation in the relativistic theory. Fluctuations inside the lightcone can propagate 
causally and tend to be positively correlated. Fluctuations in a fluid for which c$ |At| < |Ax| 
cannot have propagated from one point to the other, and are anti-correlated. This can be 
understood physically because an over density of fluid at one point in space requires an under 
density at a nearby point. 



2.2. Light Scattering by Density Fluctuations 

In Ref. [13], the cross section for the scattering of light by the zero point density fluctuations is 
computed for the case that the incident light angular frequency is large compared to the typical 
phonon frequency. The result is 
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where 9 is the scattering angle, V is the scattering volume, and rj is the mean index of refraction 
of the fluid. The u; 5 dependence of the scattering cross section can be viewed as the product of 
the u dependence of Rayleigh-Brillouin scattering and one power of Q q , and hence of u, coming 
from the spectrum of zero point fluctuations in the fluid. The factor of n represents the influence 
of the fluid on light propagation before and after the scattering process, and arises as a product 
of a factor of n in the incident flux and a factor of r/ 3 in the density of final states [13] . Because 
light travels through the fluid at speeds much greater than the sound speed, light scattering 
reveals a nearly static distribution of density fluctuations. Thus we can regard Eq. |8| as a 
probe of the fluctuations described by Eq. ©. The scattering by zero point fluctuations is 
inelastic, with the creation of a phonon. Thus, the scattering described by Eq. is strictly 
Brillouin rather than Rayleigh scattering. 

This scattering by zero point density fluctuations should be compared to the effects of thermal 
density fluctuations. The ratio of the zero point to the thermal scattering can be expressed as 
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The index of refraction, 77, and the quantity po (de/dpo) s , which involves a derivative of the fluid 
dielectric function with respect to density at constant entropy are both of order unity. Hence R 
is primarily determined by the ratio of the photon energy to the thermal energy, and the ratio 
of the speed of sound to the speed of light. 

As an example, consider the case of water at room temperature and violet light with a 
wavelength of A = 350ram. In this case, we have cs = 1480m/s and rj = 1.4 [16J. In 
addition, po (de/dpo) s = 0.79 [E]. For back scattering, cos# = —1, this leads to R « 0.005. 
Consequently, about 0.5% of the Stokes line is due to zero point motion effects. Although this 
is a small fraction, it may be detectable, and will increase at lower temperatures and shorter 
wavelengths. 



3. Squeezed States of Phonons 

Now we consider the case where the phonon field is not in the vacuum state, but rather a 
squeezed state. The squeezed states are a two complex parameter family of states, but we 
will focus attention on the case of the squeezed vacuum states |0, labeled by a single complex 
squeeze parameter 

( = re iS . (10) 



This set of states is of special interest because they are the states generated by quantum particle 
creation processes, and they can exhibit local negative energy densities. Consider the shift in 
the mean squared density fluctuations between the given state and the vacuum 

(p 2 )r = (C\p 2 \0-(o\p 2 \o), (11) 

the "renormalized" mean squared density fluctuation. The result for this quantity in a single 
mode squeezed vacuum state for a plane wave in the z-direction is 

{p 2 )n. = — 2^F smnr {sinhr — coshr cos[2(kz — cat) + 5]} . (12) 
c s V 

Note that this quantity can be either positive or negative, but its time or space average is 
positive. The suppression of the local density fluctuations in a squeezed state is analogous to 
the creation of negative energy densities for a massless, relativistic field. 



4. Boundaries 

If we introduce an impenetrable boundary into the fluid, the phonon field will satisfy Neumann 
boundary conditions 

n • VSp = (13) 

as a consequence of the impenetrability. Thus there will be a Casimir force on the boundaries 
which is analogous to the Casimir force produced by electromagnetic vacuum effects. For 
example, consider two parallel plates, which will experience an attractive force per unit area 
of 

I = (14) 
A 480 a 4 ' 1 ; 

which is smaller than the electromagnetic case for perfect plates by a factor of cs/(2c), and is 
thus quite small in any realistic situation. 

Henceforth, we consider the local effect of boundaries on mean squared density fluctuations, 
and now define {p 2 )n to be the change due to the presence of the boundary. This quantity is 
of interest both as an analog model for the effects of boundaries in quantum field theory, and 
in its own right. The shifts in density fluctuations are at least in principle observable in light 
scattering experiments. 

Our interest in the phononic analog model is inspired by the fact that the study of boundary 
effects in quantum field theory is an active area of research, and has given rise to some recent 
controversies in the literature |X8|, [T9] . One question is the nature of the physical cutoff which 
prevents singularities at the boundary. An example of the subtleties is afforded by the mean 
squared electric and magnetic fields near a dielectric interface. When the material is a perfect 
conductor, these quantities are proportional to z~ 4 , where z is the distance to the interface. 
Specifically, in Lorentz-Heaviside units their asymptotic forms are 

and 

One might expect that a realistic frequency dependent dielectric function would remove this 
singularity, but this is not the case. Instead one finds [20] that 

,;2\ \[2lO v 1 



2 > 1 ^ L , 

I 32vr z 3 v ; 



and 

where lu p is the plasma frequency of the material. Thus some physical cutoff other than dispersion 
is required. For realistic materials, it is likely to be surface roughness, but fluctuations in the 
position of the boundary can also serve as a cutoff, as we showed several years ago |21| . In a fluid, 
there is always a physical cutoff at the interatomic separation, but in a given situation other 
physical cutoffs may be dominant. For the present, we will report the results of calculations of 
(p 2 )r for several different geometries, and leave a more detailed treatment of physical cutoffs 
for later work. 

In the remainder of this section, we will quote several results for (p 2 )r in different geometries. 
The details of the calculations involved will be presented in a later paper [22J . 



4-1. One or Two Parallel Plane Boundaries 

In both of these case, the renormalized density two-point function may be constructed by the 
method of images. For the case of a single plane, one finds 

where z is the distance to the boundary. For the case of two parallel planes, the result is 

(20) 
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where a is the separation of the two planes, and z is the distance to one boundary. Note that 
both of these expression are negative everywhere. In the absence of a physical cutoff, both 
of these expressions diverge as z~ A near the boundaries, just as do the squared electric and 
magnetic fields near a perfectly reflecting plane. 

4-2. A Three- Dimensional Torus 

Here we consider a rectangular box with periodic boundary conditions in all three spatial 
directions, with periodicity lengths L±, L2 and L3. Thus the three-dimensional space has the 
topology of S 1 x S 1 x S . This is closely related to the geometry of a waveguide, where the 
fluctuations of a relativistic scalar field were discussed by Rodrigues and Svaiter [23] . As in the 
parallel plane case, an image sum method may be employed, with the result 

ipfl\ - hpo Cs nt' co-n 

2tt 2 fe n (£ 2 L 2 + m 2 L 2 + n 2 L 2 ) 2 " 1 ' 

Here the prime on the summation indices denotes that the I = m = n = term is omitted. Here 
(p 2 )r is a negative constant. 

4.3. A Wedge 

Consider two intersecting plane which are at an angle of a with respect to each other. Now 
consider a point inside of this wedge which is located at polar coordinates (r, 9), where r is the 
distance to the intersection line and 9 < a. This geometry was treated for the relativistic case 
by Candelas and Deutsch [25j, and we may transcript one of their results [their Eq. (5.39)] to 
find 

[P )R 1440^ 2 r 4 sin 4 (7r#/a) 

x {(vr - q)(tt + a) sm 2 (ir 9 /a) [(vr 2 + 11a 2 ) sin 2 (7r6>/a) - 30vr 2 ] + 45vr 4 } . (22) 



Again, this quantity is negative everywhere. 



4-4- A Cosmic String 

As is well known, the space surrounding a cosmic string is a conical space with a deficit angle 
a < 2-7T. Quantum field theory in this conical space has been discussed by many authors, 
beginning with Helliwell and Konkowski [25], and is similar to the wedge problem discussed 
above. We may follow the procedure in Ref. [25J. At a distance r from the apex, we find 

{f,2)R = - 144oX*r4 (2 " " ° )(27r + a)(Ua2 + 47f2) ' (23) 
which is also negative everywhere. 



4-5. Near the Focus of a Parabolic Mirror 

The quantization of the electromagnetic field in the presence of a parabolic mirror was discussed 
by us in Refs. [261 127] . where a geometric optics approximation was employed to find the mean 
squared fields near the focus. This treatment lead to the result that these quantities are singular 
at the focus, diverging as an inverse power of the distance a to the focus. This result holds both 
for parabolic cylinders and for parabolas of revolution, and basically arises from the interference 
term of multiply reflected rays with nearly the same optical path length. The geometry is 
illustrated in Fig. [TJ An incoming ray at an angle of 9 reflects at an angle of 9' to reach the 
point P, which is a distance a from the focus F, as illustrated. The distance from the focus to 
the mirror itself is 6/2 2> a. There can be two values of 9', denoted a and f3, corresponding to 
two reflected rays. The difference in the optical paths of these two rays is denoted by A£. 

The detailed expression for this distance Al used in Refs. [26l [27] is not quite correct, as was 
pointed out to us by Vuletic [28J. The corrected expression is 

A£ = a cos 7(cos a — cosfi + sin 2 a — sin 2 (3) + sin 7(sin a — sin (3 + sin j3 cos j3 — sin a cos a) 

(24) 

This is to be used in the expression obtained from geometric optics, 

In this expression, a = a(9), and (3 = (3(9), as will be discussed in Ref. [22]. 

A detailed discussion of the electromagnetic case will be given elsewhere. Here we are 
concerned with (p 2 )r, which is obtained from Eq. (|25p by letting c — > cs and dividing by 
2, leading to a result of the form 

(P)R= ba^ <0 - ( 26 ) 

This, and the analogous expressions for (E 2 ) and (B 2 ), which also are proportional to 1/(6 a 3 ), 
are striking in that they can be large when the focus is far from the mirror itself, 6 ^> a. 
This result is controversial, and seems to be in conflict with a general result by Fewster and 
Pfenning [29j, which implies that quantities such as (E 2 ) or (p 2 )n should be proportional to 
the inverse fourth power of the distance to the mirror, which is to say oc 6~ 4 in this case. On 
the other hand, there is a simple physical argument to the contrary in this case, which we 
find compelling: the interference term between multiply reflected rays is slowly oscillating when 
A£ oc a is small, and should give a contribution proportional to an inverse power of a, as in 
Eq. (j26l) . In any case, the study of the phononic case provides an additional theoretical, and 
potentially experimental, probe to better understand this issue. 



Figure 1. The geometry of rays reflecting from a parabolic mirror is illustrated. An incoming 
ray at an angle of reflects at an angle of 9' to reach the point P, which is a distance a from 
the focus F, and at an angle of 7. 

5. Summary 

In this paper, we have considered a classical or quantum fluid with a linear dispersion relation 
as a analog model for quantum fluctuations, and the effects of boundaries on these fluctuations. 
We have argued that the local density fluctuations in a fluid in the phonon ground state are 
potentially observable in light scattering experiments. This is of interest in its own right. We 
have reported the calculation of the change in these fluctuations due to squeezed states of 
phonons and due to the presence of boundaries. Squeezed states can produce either positive or 
negative values for (p 2 )r, in a way which is analogous to the effects of squeezed states on the 
energy density of a relativistic quantum field. However, all of the boundary examples which 
we consider result in a negative (p 2 )r- If we compare to the already negative result in Eq. (J7|), 
this seems to represent an increase in magnitude of an already negative quantity. If one were 
to scatter light from the fluid with a boundary, this implies a change in the scattering cross 
section, compared to the result for a boundary- free fluid, Eq. ([8]), at least for short wavelengths. 
In general, the calculation of the modified cross section in the presence of a boundary is a complex 
task. However, if the incident wavelength is short compared to the scale over which {p 2 )u varies, 
then we might the scattering cross section to increase in the presence of a boundary. A more 
detailed analysis of this effect is a subject for future study. It will also be of interest to consider 
the effects on non-linearity in the phonon dispersion relation on the questions considered here. 
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